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PREFACE 


Atmospheric  models  are  frequently  generated  for  specific  regions  of  the 
atmosphere  in  accordance  with  a  particular  set  of  observations.  For  example, 
for  the  region  20  to  90  km  or  for  the  region  120  to  100  km  altitude.  These 
models  are  not  necessarily  extended  to  sea  level,  neither  are  they  designed 
to  be  continuous  with  some  other  existing  model  at  either  end  of  the  included 
range.  In  some  instances  it  becomes  desirable  to  exactly  connect  two  such 
models  with  a  transition  model  atmosphere  which  for  the  examples  cited  would 
extend  between  90  and  120  km.  This  report  discusses  a  method  suitable  for 
generating  such  transition  models. 
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ABSTRACT 


Any  set  of  temperature , pressure  and  density  values  which  are  realistic 
for  one  altitude  can  be  exactly  connected  to  another  set  of  temperature, 
pressure  and  density  values  at  a  second  altitude,  with  a  model  atmosphere 
defined  by  the  appropriate  linearly  segmented  two-layer  temperature-altitude 
profile  provided  one  of  these  layers  is  isothermal  and  the  second  is  a  nonzero 
constant  gradient  of  temperature  with  respect  to  height. 

The  altitude  comprising  the  intersection  of  the  two  segments  of  the 
required  temperature-altitude  profile  is  determined  mathematically .  The  method 
for  generating  two-layer  models  is  applied  to  the  generation  of  three-layer 
and  four-layer  models. 
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SECTION  I 


REVIEW  AND  EXTENSION  OF  THE  METHOD  FOR 
DEVELOPING  TWO-LAYER  TRANSITION  MODELS 


l.  DEFINITION  OF  THE  GENERAL  TWO-LAYER  TRANSITION  MODELS  AND  THE 

SPECIFICATION  OF  THE  PARTICULAR  CLASS  TO  BE  EXAMINED 

Minzner  (1966,  pp.  5  and  6)  has  shown  that,  under  a  particular  though  not 
highly  restrictive  set  of  conditions,  a  two-layer  transition  model  can  be  made 
to  fit  exactly  between  two  sets  of  fixed  boundary  conditions.  The  first  con¬ 
sists  of  a  set  of  realistic  values  HQ,  p0  and  70,  where  H0  is  the  geopotential 
corresponding  to  one  arbitrary  geometric  altitude,  p0  is  a  corresponding  rea¬ 
listic  value  of  density,  and  Tc  is  a  realistic  value  of  the  temperature  para¬ 
meter  at  H0.  The  second  consists  of  a  set  of  realistic  values  1^,  pn  and  Tn, 
where  Hn  is  the  geopotencial  corresponding  to  another  arbitrary  geometric 
altitude,  Pn  is  a  corresponding  realistic  value  of  density,  and  Tn  is  a  cor¬ 
responding  realistic  value  of  the  temperature  parameter.  The  subscript  im¬ 
mediately  to  the  right  of  the  symbols  H,  p,  or  T  is  associated  with  the  par¬ 
ticular  geopotential  level  to  which  these  three  quantities  are  related.  In 
this  discussion  of  transition  models,  the  words  altitude  or  height  will  be 
intended  to  imply  the  equivalent  geopotential  without  such  a  statement  neces¬ 
sarily  being  made. 

The  restrictive  conditions  defining  the  particular  class  of  two-layer 
models  to  be  considered  are  as  follows: 

(1)  One  of  the  two  layers  is  characterized  by  a  non-zero  constant 
gradient  of  the  temperature  parameter  relative  to  the  geopotential  such  that 
this  non-zero-gradient  layer  extends  from  H0  toward  Hn  for  an  unspecified 
altitude  interval  to  the  level  H^. 

(2)  The  second  layer  extending  from  to  1^,  is  characterized 

by  a  constant  zero  gradient  of  the  temperature  parameter  relative  to  the  geo¬ 
potential.  The  symbol  Hxo  designates  the  altitude  of  the  interface  between 
the  isothermal  layer  and  the  non- isothermal  layer  extending  up  to  H0. 

The  equations  and  computat ional  procedures  for  generating  the  class  of 
transition  models  as  discussed  in  this  report  are  developed  on  the  basis  of 
two  further  assumptions:  (1)  Geopotential  Hn  represents  a  lower  altitude 
level  than  that  represented  by  Hc  such  that  the  isothermal  layer  is  the  lower 
of  the  two  layers  of  the  transition  model,  and  (2)  the  value  of  the  tempera¬ 
ture  decreases  with  decreasing  altitude  from  H0  to  the  isothermal  layer.  This 
class  of  model  applies  to  situations  similar  to  that  portion  of  the  W.  S. 
Standard  Atmosphere  1962  between  altitudes  of  117  to  79  !  ilometers. 

The  temperature  parameter  employed  in  the  transition-model  development  is 
the  molecular  scale  temperature  which  combines  two  variables,  the  kinetic 
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temperature  and  the  mean  molecular  weight  of  the  atmosphere ,  into  a  single 
variable  in  accordance  with  the  following  definition: 

The  molecular  scale  temperature,  at  any  level  II,  is  equal  to  the 
ratio  of  the  kinetic  temperature  to  the  mean  molecular  weight,  at  that  same 
level,  times  Ms  the  sea-levol  value  of  the  molecular  weight. 

For  simplicity,  the  notation  commonly  used  for  molecular  scale  tempera¬ 
ture  will  be  avoided  in  favor  of  T  in  this  discussion  of  transition  models, 
and  the  word  temperature  will  imply  molecular  scale  temperature. 

2.  DENSITY  AND  TEKPERATl'RE  RELATIONSHIPS  FOR  THE  SPECIFIED  TYPE  OF 

TWO- LAYER  MODEL 

In  the  type  of  model  hypothesized,  where  an  isothermal  layer  at  tempera¬ 
ture  Tn  extends  upward  from  Hn  to  Hxo  and  a  layer  with  a  constant  non-zero 
temperature-altitude  gradient  extends  upward  from  Hxo  to  !!0,  the  density  0X 
at  Hx  may  be  expressed  in  terms  of  on  and  Tn  by  the  following  equations: 


-(H  -  H  )  Q' 

xo  n  (.  , , , 

o  =  q  exp  - r  (i) 

xo  n  r  T  J 

■-  n  J 

where 

C  M 

Q  *  =  0.034J  6319470K/ui'  (2) 

K 

and 

i:  is  the  geopotential  transformation  coefficient  equal  to 

9.80b65  sec'2  (m’)“l 

M.  is  the  sea-level  value  of  the  molecular  weight  of  air, 

S  28.9b44  kilograms  per  kilomole  (kg/kmol) 

R  is  the  universal  gas  constant  8.31432  x  10^  joules  (°K)-* 
(kmol)~l 


In  such  a  model  the  density  ox  at  Ilx  may  also  he  expressed  relative  to 
Tn  and  the  upper-level  boundary  condition  values,  pQ,  T0,  and  H0,  which  through 
Hx  imply  the  value  of  the  non-zero  gradient  of  T  between  Hx  and  H0.  Thus , 
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The  value  of  Hxo  is  found  by  the  simultaneous  solution  of  Equations  (1) 
and  (3)  as  expressed  by; 


The  temperature-altitude  gradient  between  Hxo  and  HQ  specified  by  L^ 
is  given  by 


L 


xo 


(5) 


Here,  the  symbol  for  the  gradient  has  two  subscripts,  the  first  designating 
the  unknown  lov"r  end  of  the  associated  layer  and  the  second  designating  the 
known  upper  end  of  the  associated  altitude  interval. 


3.  ANALYTICAL  EXPRESSION  FOR  THE  GENERATION  OF  THE  TWO-LAYER  MODELS 

Because  It  was  not  at  first  suspected  that  an  analytical  solution  for 
Hxo  in  Equation  (4)  could  be  found,  numerical  methods  involving  digital  com¬ 
puters  were  used  to  determine  the  value  of  Hxo  for  a  number  of  transition 
atmospheres  (Minzner,  1966,  pp.  7  to  24).  It  may  now  be  shown,  however,  that 
an  analytical  solution  for  Hxo  in  Equation  (4)  does  exist. 

Dividing  both  sides  of  Equation  (4)  by  0o  and  taking  the  natural  logari- 
thum  of  both  sides  of  the  resulting  expression  leads  to 
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The  separation  of  the  numerators  of  the  two  fractions  containing  Hx,  so  as  to 
permit  the  extraction  of  Hxo,  followed  by  a  rearrangement  of  the  terms  leads 
to  tiie  expression 
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Factoring  Hxo  out  of  the  left-hand  side  oi  ttu  equation  .nut  d  i  v  id  iin;  both 
sides  of  that  equation  bv  the  remainder  of  the  left-hand  side  yields 


An  evaluation  of  this  expression  using  the  boundary  conditions  of  model 
Ar  (Minzner  1966,  Table  2,  p,  13)  yielded  Hxo  =  101.9646  km.  This  value  agrees, 
to  within  one  part  in  the  sixth  significant  figure,  with  the  value  obtained  by 
the  numerical  methods  (Minzner,  1966,  Table  3,  p.  14),  and  thereby  simultaneously 
substantiates  the  validity  of  the  numerical  method  as  well  as  of  the  analytical 
solution  of  !lxo  given  by  Equation  (8). 

The  subscript  o  to  the  right  of  x  in  Hxo  is  used  to  indicate  that  this 
particular  value  of  Hx  is  calculated  with  the  upper  boundary  conditions  at 
level  H0 .  If,  as  will  bo  seen  in  the  generation  of  a  three- layer  model,  Hx 
is  calculated  from  some  other  upper-level  boundary  condition,  as  at  the 

notation  associated  with  x  will  be  m  as  in  H^. 

The  value  Hxo  1  rom  Equation  (8),  and  the  value  Lxo  from  Equation  (5), 
along  with  the  upper  and  Lower  set  of  boundary  conditions  completely  define 
the  only  two- layer  transition  model,  of  the  type  under  consideration  which 
satisfies  the  boundary  conditions.  There  are  multilayer  models  (three  or 
more  layers),  however,  which  satisfy  the  same  upper  and  lower  set  of  boundary 
conditions  and  the  unique  two-laver  model  serves  as  a  basis  for  the  generation 
of  these  multilayer  models.  Computer  Programs  for  the  generation  of  two-layer 
models  arc  given  in  Appendix  A. 
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SECTION  II 


GENERATION  OF  THREE -LAYER  MOTELS 


The  method  for  the  generation  of  a  two-layer  model  between  upper  and 
lower  boundary  conditions  can  be  extended  to  the  generation  of  an  infinite 
number  of  three-layer  models.  The  method  consists  essentially  of  arbitrarily 
selecting  an  altitude  interval  for  the  uppermost  layer,  computing  the  conditions 
at  the  base  of  that  layer  from  the  initial  upper-boundary  conditions,  and  using 
these  new  values  as  the  basis  for  the  computation  of  the  related  two- layer 
model.  The  detailed  steps  to  be  used  are  as  follows: 

(1)  Compute  the  value  of  Hx0  and  I^q  for  that  two-layer  model 
which  fits  the  initial  boundary  conditions. 

(2)  Select  a  suitable  value  as  the  base  of  the  uppermost  or 

first  layer  of  the  desired  three-layer  model.  A  suitable  value  meets  the 
condition  Hxo  <  <  H0,  An  infinite  number  of  possible  values  of  Hj  exist 

between  Hxo  and  H0. 

(3)  Select  a  suitable  value  of  for  the  layer  HQ  to  Hj.  It  may 
be  demonstrated  that  if  the  gradients  of  the  several  layers  are  to  decrease 
monitonically  from  H0  to  Hn,  the  value  of  must  meet  the  condition  Lxq<  L^. 

Furthermore,  it  will  be  seen  that  Lj  must  meet  the  condition 

L1  <  Imiax  "here  L,,^  is  an  upper  limit  associated  withTlminat  Hj  as  in  Figure  1. 

The  determination  as  to  whether  Li  <  is  based  upon  the  result  of  a  test 

following  a  later  step  in  the  proceedure.  At  this  point,  however,  a  trial 

value  of  is  chosen.  An  infinite  number  of  possible  values  of  meet  the 

necessary  condition  <  L,  <  L 

7  xo  i  max 

(4)  Compute  at  the  selected  value  from  the  selected  value  L^ 
using  the  relationship 


T1  ■  To  -  Ll  (Ho'Hl>  (9) 


(5)  Compute  the  value  of  p^f  the  density  at  H  ,  using  the  rela¬ 
tionship 


\  "  Po 


-ii  + 


<Ho  -  Hi>  q  I 
To  -  T1  j 


'10) 
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The  set  ot  values  llj,  Tp  and  ,  j  serve  as  a  new  upper- boundary  condition  which 
with  1^,  rn,  and  P  provide  the  means  for  generating  a  two-layer  model  between 

111  and  H  ,  where  these  two  layers  are  designated  layers  2  and  3  of  the  three- 

layer  model.  The  Interface  hetween  layers  2  and  3  Is  designated  Hxj. 

<bi  Determine  Hx^  using  Equation  (8)  with  suitable  subscript  modi¬ 
fications.  i.e.,  subscript  zero  is  changed  to  one.  If  has  been  properly 
selected  In  step  3,  llxi  will  be  found  to  have  a  value  such  that  "n  '  Hxl  <  “xo* 

and  this  condition  serves  as  the  test  for  suitability  of  I,j, 

t'l  Compute  the  gradient  for  the  layer  between  llx,  and  by  using 
Equation  t"»3  with  the  subscript  zero  Increased  to  subscript  1. 

Vhe  ,  .'millions  for  the  test  of  suitability  of  hj  based  on  the  value  of 
u  have  been  developed  through  the  following  argument: 

It  is  apparent  that  if  is  selected  to  he  equal  to  l^p,  11^  is  not 

real l v  the  base  of  a  new  layer,  but  is  a  level  within  the  layer  Hx0  to  !!„. 

Computations  based  on  lllt  Tj,  and  .  i  would  then  yLeld  llxj  -  Hxo ,  and  one  would 
have  determined  only  a  redundant  case  of  the  original  two-layer  model.  As  the 
value  of  l.t  is  allowed  to  Increase  from  l^0,  for  a  fixed  value  of  Hi,  the  value 
of  H  .  progressively  decreases  from  Hxo  until,  in  the  limit  for  Li  “  l^nax* 

value  of  :lxl  becomes  identically  Hn  in  which  case  the  isothermal  layer  has 

vanished  as  indicated  In  Figure  1.  For  values  of  Li  >  L  ,  the  computed 

value'  of  'l  become  less  than  11  ,  indicating  an  impossible  situation. 

xl  n 

It  is  also  possible  to  select  values  of  Iq  such  that  I^in  <  <  ^to" 

As  Lj  Is  decreased  from  the  value  of  Hxi  increases  trom  llxo  until 

for  4,1,,,  llxi  *  Hj.  For  this  case,  layer  2  has  the  impossible  conditions  of 
zero  thickness  and  an  infinite  temperature  gradient.  All  values  of  Li  < 
are  Impossible  however,  If  the  model  is  constrained  to  one  having  a  monl- 
tonical ly  decreasing  set  of  values  of  L  for  successive  layers.  It  follows 
therefore,  that  the  only  suitable  values  of  Lj^  are  those  which  lead  to  values 
of  Hxl  such  that  <  Hxl  <  Hxo.  An  infinite  number  of  such  values  exist. 

It  values  of  Lj  are  constrained  to  be  integral  multiples  of  some 
fixed  values  J-i ,  and  if  they  are  constrained  to  be  within  the  range  to 

I.  .  the  smallest  suitable  values  of  L,  may  be  expressed  as 
max  * 


L 


la 


h  -U1 


(ID 


where 


satisfies  the  inequality 


<  »i  <  >  + 


(12) 
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GEOPOTENTIAL 


Figure  1.  Schematic  temperature-altitude  diagram  for  a  two-layer 
model  atmosphere,  for  an  acceptable  related  three-layer 
model,  and  for  a  limiting  case  of  the  three-layer  models. 
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The  series  of  acceptable 

11  ,  are  then 

x ! 

values 

of  I.j 

and  the  related 

series  of 

Lu  •  <Ji  ♦ 

1  -  1) 

U1  ’ 

H  i 

xla 

(13a) 

ho  ■  <h  * 

2-1) 

M'l  • 

H 

xlb 

U3b) 

he  ’  «l  f 

3  -  t) 

**1  ’ 

!!  t 

XIC 

(13c) 

etc.,  until  a  value  Ljg  is  reached  for  which  is  less  than  Hn.  This  sit¬ 

uation  is  shown  schematically  in  Figure  2. 
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0SCAI92*  90S 


GEOPOTENTIAL 


Figure  2.  A  schematic  temperature-altitude  diagram  showing  four  possible 
three-layer  models  and  one  impossible  three-layer  model. 
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SECTION  III 


GENERATION  OF  MULTILAYER  MODELS 


The  above  procedure  nay  be  extended  to  the  generation  of  four-layer 
models  and  on  to  the  generation  of  n-layer  models.  The  procedure  consists  of 
generating  a  three-layer  model  as  indicated  above,  and  then  taking  the  lower 
two  layers  of  this  model,  follow  the  procedure  of  steps  2  through  7  of  Section 
II,  but  with  subscript  l  replaced  by  subscript  2,  and  subscript  o  replaced  by 
subscript  1.  In  a  similar  manner,  the  procedure  may  be  extended  to  a  model  of 
any  number  of  layers. 

The  altitudes  specifying  the  boundaries  of  the  four  layers  of  each  of  a 
set  of  several  possible  four- layer  models  for  a  specified  set  of  boundary 
conditions  are  given  in  Table  1. 
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APPENDIX  A 


Program  1 


The  program  for  defining  a  two- layer  transition  model  atmosphere  by  the 
meched  of  successive  approximation  and  from  the  definition  computing  the 
various  properties  as  indicated  in  Tables  3  through  10  of  the  final  report 
or.  Contract  No.  NASV-1225  (Minsner,  1966). 


Preceding  Page  Blank 


PROCRAM  l 


source  DFr*  FOR  THF  c  DMpi  iTATI  ON  OF  TWO-LAYfR  transition  MODEL 

ATMOSPHFRF  A  RFTwFFN  tIPPFR  t.AYFR  AND  LOWFR  (AyFR  ROUNDARy  CONDITIONS 
WHFW  TmF  cn|  CONDITIONS  REGARDING  TFmPFR  A  Tuff F— AL  T  I  TUDF 

GRADIENT  PPFVA  f  (  S 

T ,  IT  MAS  A  O I FF  FRF  NT  CONSTANT  VALUE  FOR  FACh  iAyFR. 

?,  it  has  a  positive  valuf  for  the  uppermost  layer. 

3.  IT  HAS  A  ZFRO  VALUE  FOR  THF  LOWEST  LAYER. 

THF  METHOD  INVOLVES  THE  SI  MULTAN  I OUS  NIJMFR I  CAL  SOLUTION  OF 
TWO  EQUATIONS  USING  THF  MFTHOO  OF  SUCCESSIVE  APPROXIMATIONS 
as  0 1 SCHSSFD  BY  MINZNFR  IN  GCA  TECH  RPT.  N0.66-14-N 
J<  l\iF  1QAS,  RFVISFH  M4RCH  1966 
’•»  RC40  1^1  .0FG.HB.H»R.PMP«HA.HP4,PHA,FM0N.T1,T2 

READ  IN  INITIAL  CONDITIONS  FOR  MODEL 
T 1  AND  T?  ARE  THF  TITLE  FOR  THE  MODEL 
•  -  *  r^pMftr  (4-».F4#D*F10.7,F11.A.F9.4.F9,S«F10.1,1XA4.?AS> 

QS  An  IIS.FL.HRP 
T  1  s  FORM4T<F/».  1  .  1XFS.01 

FPSJ=.1F-0R 
X  •  1 

PUNCH  11? 

n?  FORMATI3SHPROPFRTIFS  OF  TRANSITION  ATMOSPHERE  » 

PUNCH  106.T1 .T? .DFG.FMON 

106  FORMAT(?AS.2xlH( .A3.1HI »ZXA4,4X17HT0  JACCHJA  MODELS/) 

SOLVE  THE  TWO  EQUATIONS  SIMULTANEOUSLY  USING 
THF  METHOD  OF  SUCCFSSlVF  APPROXIMATIONS 

7  PFTA=(HA-X 1 /(HPA-HPP) 

Y  -  f  1  ,+RF  T  A  )  *LOC>  (HP  A  /HPBJ+LOG  (  PH  A  1 
xni  D*x 

X  *-  I  Y-LOG  (  PHB  )  )  *HPR  /LOG  (  E  )  +HR 
IFfSFNSF  SWITCH  H  4.3 

SENSE  SWITCH  1  ON  WILL  LET  THE  OPERATOR 
WATCH  THF  D I  RFC  TON  THE  X  VALUE  IS 
GOING  IN. 

i  POINT  107.x 
inx  FORMAT ( F 14. fl  ) 

4  IF ( ARSfxOLO-X l-FPST  11  .  1  ,  ? 

1  CONTINUF 

Y  * FXP ( Y ) 

PFTA=1./RFTA 

XOT=X+0. 

PUNCH  10?»DFG.FM0N.X  .Y 

107  FORMAT f A3 , 1 XA4, ?*3HX  «F 1 0 . S * ? X ?H<M4X 7HDFN  X  =F 14.8 .2X4HG/M3  ) 

PUNCH  1 1 3.HA.HPA.HPP 

113  FORMAT ( 19HGFOPOT  SCAlF  HT  AT  F9.4.1X4HKM  =F9.5.1lH»  AT  X  KM  =. 

1  F  1  0 . 7  ) 

PUNCH  1D9.PETA 

1PQ  FORMAT ( 3SHGRAD! FNT  OF  GEOP  SCALE  HT  ABOVF  X  *E14.8) 

Punch  H4.FL.HRP.hr 

114  FORMAT ( ?3HGRAD1 FNT  OF  MOL.  TEMP.  .F4.1.10H  FOR  BTWN  .E5.0.4H  TO  »F4.0) 
Punch  IIO.mB 

11'*  FORMAT  {  XAHORAO!  FNT  OF  MOL.  TEMP.  ZERO  “OR  H  RTWN.F4.0.6H  AND  X) 

punch  107 

*07  FORMAT* 3X 1 ’HALT 1TUDF  KM10X1SHTEMP  DEG  KEL V I N4X 6HMOL  WTZX. 
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117HDFNMTY  G/M34X  1 1  HSCALF  HT  KM) 

PUNCH  108 

1  '">«  FORMAT!  lX6HGFOMFT4X6HGFOPOT7X4HMOL«3X7H<iNFT IC28X6HGEOMET1X. 
RHOX  =  V 

RPAD  104,DFG,R  ,GPHI 
1  04  FORMAT*  A3»F9.0»  FP.8  ) 

H  =  XPT 

Z  =  R*H / ( GPH I *R/9 .80665— H I 

C  Z  IS  THE  GEOMETRIC  ALTITUDE 

RHOH*PHA*(HPA/(MPA+RETA*<H-HA) ))**!( 1 .+RFTA ) /BFTA) 

C  RHOH  IS  THE  DENSITY 

HH=  ( HPR+ BETA*  ( H-X  )  ) *9.80668 /GPHI* f  !R+7)/R)*»? 

C  HH  IS  THE  SCALE  HFTGHT 

TM=34.163?*(HPB+BFTA*(H-X) ) 

T  =  TM*( 1.-2.37S81 lE-03*(Z-90. ) ) 

C  T  IS  THE  TEMPERATURE 

HPH=HPR*RFTA* ( H-X ) 

SM=28.O644-.06881 33*fZ-90.  ) 

C  SM  IS  THF  MOLFCULAR  WEIGHT 

TM=TM+,nOS 
T  =  T  +  .008 
SM=SM+.OOS 
HH=HH* . 0008 
HPH=HPH+.0005 

C  ROUND  THE  VALUES 

PUNCH  108*Z*H,TM*T*  SM . RHOH  *HH  »HPH 
C  PUNCH  DATA  CARD 

IRS  FORMAT f F8.3t?XF8.7,3XF7.?.3XF7.2,3XF6.2.?XF12.5 .2XF7 .3 .2XF7 . 3  I 
16HGFOPOT/ ) 

C  METHOD  FOR  FINDING  NFXT  EVEN  KM. 

XNFW=DRH  *  X  I 
XNFW»XNFW/2. 

XNEW=DRH ( XNEW ) 

XNFW=I XNFW*2. >  +  2. 

XP=X*10. 

XP=ORH(  XP ) /10. 

IF! XNFW-XP-7. )  8  »9 ♦ 8 
o  H  =  XP 
GO  TO  7 
8  H=XNEW 

7  I F (H—XPT )  5  »5  ♦?  3 
?3  Z=R*H/(GPHI*R/9 .80668-H) 

RHOHspHA* I  HP A/ ( HP A+RET  A*  t  H— HA ) ) ) ** C  t 1 .+BETA ) /BETA ) 

HH= (HPB+RFTA* ( H-X ) ) *9 . 8066 5/GPHI * C (R+Z )/R)*#2 
TM=34.163?*(HPB+BFTA*(H-X) ) 

T=TM*< 1 .-?.375811F-03*CZ-90. ) ) 

HPH=HPB+8FTA* ( H-X ) 

SM*28 • 96 44 -.06 8 8 133*(Z-90. ) 

FLM=34.1632*BETA 

C  FLM  IS  THE  GRADIENT 

TMaTM+.OOS 
T*T+.008 
SM*SM*.008 
HH*HH+.0008 
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HPWsHPH+.0/'05 
IF  f  7~!70.  U0.6«6 

40  PUNCH  106»7*H»Tm«T  «  SM  .RHOH  *HH  »HPH 
*1  H*H  +  7* 

on  TO  7 

4  7rO*X/rCPHf  *R/<?  .80665-X  ) 

7 sDPH? 7  J /2. 

z=orh<z  »#?.+?. 

17  h-*R*7/(R+Z)*GPH!/9. 80665 

RH0H=PHA* (HPA/ ( HP A ♦BE  TA*  ( H-HA ) )  )  **  ( ( l.+BETA  J /BETA  I 
<,M  =  ?8,Q644-.06  6813  3*<  Z-90. 1 

HH=(HPR+8FTA#IH-X 1 >#9.80fc65/,GPH!*{ {R+Z)/R>#*2 
Tv=34. 163  7*  <HPB+BFTA* (H-X  > 1 
T=TM*( l .-?,37A81ir-03*I7-90. ) > 

ELw-74. 163?*f)FTA 

HPH=HpR+BFTA#(H-X I 

TM=TM+.005 

T=  T  +  .005 

6m=6M+.005 

HH=HH+.0G0S 

HPH=HPH+.00C5 

PUNCH  105  » 7  .H.TM.T ,SM,RHOH*HH.HPH 
IF  (7—170.  HO, 11  ,11 
10  7=7+7. 

CO  TO  17 

1  1  rlsHR 
X  =  HR 
BFTA=0. 

10  7  sp#H/(  C,PH  r  *R/9 .8066  5-H  I 

HH=(HP«  +  RFTA*( H-X ) J  *9 . 80665 /G PH I  * (  (R+2 ) /R  >**? 

TMs3A. 1637* (HPR+BFT4* (H-X I ) 
j F  f 7-90. ) 26.76 .77 

76  T=TM 
FMs?8.96 
GO  TO  78 

77  T=TM*( 1 .-7.37581 1E-03*(Z -90. ) > 

?3.  Oft  a*-.  06  881  3-»*(7-90.  ) 

70  dhoh  =  PHP*FXP I  —  I H— HR ) /HPR ) 
mDH«HPB+8FTA*{H-X ) 

XPPlRK=DRHr  XPT  ) 

Tm=TH+.005 
T  =  T  +  »005 
4m=4m+. 005 
HH=HH+.0005 
hph=HPH+«0005 
IF  IH-XPRIMHA,!  A. 15 
I*  punch  io6»z»h,tm,t.4m,rhoh.hh,hph 
IF?H-HR) 16.41,16 

41  I F ( SFNSF  SWITCH  7)16,17 

SENSE  SWITCH  2  ON*  INCREASE  THE  HEIGHT  BY  2. 
SENSE  SWITCH  2  OFF,  INCREASE  THE  HEIGHT  BY  1. 

1 7  H=H+1 • 

GO  TO  18 
16  H=H+7. 
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GO  TO  18 

1  8  2 sR*X / ( GPHI *R /9 • 80665-X  J 
Z=ORHCZl/2. 

2*DRH<Zl*2.+2. 

7 1  H=R*Z/fR+Z)*GPHI/9. 80665 
I F (H-HB) 25 .25  *2? 

77  RHOH»PHB*FXP< -(H-HB) /HPB  1 

HH*(HPR+BFTA*IH-X1 1*9 *80665 /GPHI *t {R+Z1/R)**2 
TM=34, 1632* ( HPB+BET A* ( H-X 1 1 
IFfZ-90. 129*29. BO 

20  T  =  TM 

$M*28.96 
GO  TO  31 

30  T  =  TM*( l.-2.375811F-O3*«Z-90.> 1 
Fm=28.9644— *0688 1 33* ( Z  —90* 1 

n  HPH*HPP*8FTA* (H-X 1 
TMaTM+,005 
T=T+*005 
6M*SM+*005 
HH=HH+«  0005 
HPH=HPH+*  0005 
IFfH-XPTl 19*20.20 

1 Q  PUNCH  105»Z»H»TM*T  *  5M  »  RHOH  *HH  *HPH 

26  1*7*7. 

GO  TO  21 

20  PUNCH  111 .FLM.HA 

111  FORMAT ( 23HGRADIFNT  OF  MOL.  TEMP.  F10.5 .2X19HDEG/ICM  FOR  H  BTWN  X. 

14H  AN0.FQ.4//1 


PROGRAM  2 


The  program  for  defining  a  two- layer  model  by  the  evaluation  of  an 
algebraic  function,  which  replaces  successive  approximation  method. 
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rs  r\  pi  n  r\ 
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PROGRAM  2 


SOURCF  DECK  FOR  TWO-LAYER  TRANSITION  MODEL  ATMOSPHERE 
^«Ff^°^PfR"LFVEL  AND  LOWER-LEVEL  BOUNDARY  CONDITIONS 

f°LLOWING  CONDITIONS  REGARDING  TEMPERA TURE-ALT I TUDE 
GRADIENT  PREVAILS 

I.  IT  HAS  A  DIFFERENT  CONSTANT  VALUE  FOR  EACH  LAYER. 

?•  IT  HAS  A  POSIT  I VF  VALUE  FOR  THE  UPPERMOST  LAYER* 

’*  IT  HAS  A  ZERO  VALUE  FOR  THE  LOWER  LAYER. 

*1?  F^A^F12-1*£xe^^.lXF10.3*2xFin.l.2XE12.6*lXF10.3» 

51  READ  100  *HN .RHON.TN.HO .RHOO  *  TO 
BEGIN  TRACE 
Q*3.4l6319«>E-02 

HXn. ( LOGF I RHON /RHOO ) +HN*Q/TN-  ( 1 . +H0*Q/ 1 TO-TN I  I *LOGF I TO/TN  >  > 
HXO»HXO/(Q/Tn-IO/(TO-TN) >  *LOGF ( TO/TN ) I 

GLXOs( TO-TN J/IHO-HXO) 

RHOXO*RHON*ExPF I  -o*  < HXO-HN  J /TN ) 

120  FORMAT ! 2 1X33HDEF IN  I T ION  OF  THE  TWO-LAYER  MODEL  I 
111  FORMAT ( 1H  ) 

11?  FORMAT ( 10X5HLEVEL  *  5X 1HH *  1 1X3HRH0 *  1 IX 1HT  *  1 IX 1HL ) 

113  FORMAT! 18X43HMETERS  K6/M3  DEG  K  DEG /M ) 

114  FORMAT! 10XI3»2XF10.1.2XE12.6*1XF10.3.?XE12.6) 

1  T<5  nf!^uT!  !?XI3,?XF10*1  *^xF^2.6.1XF10.3»3X  1  lH.OOOOOOE+OO ) 

PI  fNCM  111 
PUNCH  HI 
PUNCH  120 
PUNCH  HI 
PUNCH  112 
PUNCH  113 
PUNCH  111 

lfvfl=o 

PUNCH  114.LEVEL.H0.RH00.T0 
LFVEL*LEVEL+1 

PUNCH  1I4.LEVEL.HXO.RHOXO.TN.GLXO 
LEVEL*LEVEL+1 

PUNCH  115 .LEVEL. HN .RHON »TN 

END  TRACE 

END 


